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Abstract. This paper proves a theorem about bounding orbits of a time 
dependent dynamical system. The maps that are involved are examples in 
convex dynamics, by which we mean the dynamics of piecewise isometries 
where the pieces are convex. The theorem came to the attention of the authors 
in connection with the problem of digital halftoning. Digital halftoning is a 
family of printing technologies for getting full color images from only a few 
different colors deposited at dots all of the same size. The simplest version 
consist in obtaining grey scale images from only black and white dots. A 
corollary of the theorem is that for error diffusion, one of the methods of 
digital halftoning, averages of colors of the printed dots converge to averages 
of the colors taken from the same dots of the actual images. Digital printing 
is a special case of a much wider class of scheduling problems to which the 
theorem applies. Convex dynamics has roots in classical areas of mathematics 
such as symbolic dynamics, Diophantine approximation, and the theory of 
uniform distributions. 



1. Introduction 

Theorem 1.1. Let P be a polytope in , and for each j in P let (j)^ denote the 
map from to defined by 



(1.1) (f>^{x) =x + {j-v{x)), 

where v{x) is the closest vertex of P to x in Euclidean norm with some tie-breaking 
rule. Then for any compact set K there exists a bounded convex domain Q contain- 
ing K and invariant in the sense that for all j in P, (j)~fQ <Z Q . 

We encountered Theorem ll.ll in connection with a problem in digital color print- 
ing 0. It turns out to have wide applicability (see, e.g., ^7]) and considerable 
mathematical depth. 

The general problem that this theorem addresses concerns approximating a 
bounded one-sided sequence 7(fc) of arbitrary values in by a sequence of el- 
ements V{k) chosen from some finite set. In a wide range of applications, such 
as digital printing or scheduling, a reasonable sense of well-approximation is that 
there be a uniform bound (small if possible) on some norm of the cumulative error 
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vectors 

n 

(1.2) s{n)^Y.i^ik)-V{k)), 

in which case the average error wiU go to zero. In apphcations it is assumed that 
the V{k)^s are taken as the vertices of a polytope P and the 7(fc)'s belong to this 
polytope. In the case when the norm 1 1 • 1 1 is the sup norm and P the standard 
simplex, the problem of finding the best universal bounds on sup„ ||e(?T-)|| is often 
referred to as the Chairman Assignment Problem. This problem was posed by 
Niederreiter in \2'2\ and considered in ,,20ul9t y33u34) : we say Generalized Chairman 
Assignment Problem when the norm is not necessarily the sup norm. Using the 
methods of JJj to find sequences that solve the Generalized Chairman Assignment 
Problem turns out not to be effective as it requires looking ahead with unbounded 
waiting time, at least for the sup norm. Waiting for future input data before making 
output decisions is not practical in many cases and sometimes even not acceptable 
as the scheduling algorithms may need, not only to be causal (i.e., depend only 
on the past and the present and not also the future) but also to be on the fly 
(i.e., provide outputs essentially as soon as input data are known). Thus, if one 
is not looking for the best bound but only for a uniform one, or if one needs to 
perform continual optimization, then a surprisingly good approach is to write e{k) 
recursively as 

(1.3) e{k + 1) = e{k) + 7(fc + 1) - v{e{k) + ^{k + 1)) , 

and use the greedy algorithm to specify v{e{k) + 7(fc + 1)) which means here that 
v{£(k) + 7(fc + 1)) is the vertex which minimizes the norm of e{k + 1) with some 
tie breaking rule: i.e., the vertex closest to e{k) + 7(fc + 1) with some tie breaking 
rule. We shall be interested in the Euclidean norm; and in what follows, when 
we say "the" greedy algorithm, we mean the greedy algorithm using the Euclidean 
norm with any tie breaking rule for points equidistant from two or more vertices. 
Later it will be convenient to express the greedy algorithm in terms of the notion of 
Voronoi' regions. A Voronoi region of a vertex is the closure of the set of points that 
are closer to the vertex than to any other vertex (see Definition 12. 1(1 . Alternatively 
expressed, the greedy algorithm chooses the vertex in the Voronoi region containing 
e{k) + 7(/c + 1) with some decision rule when e{k) + ^{k + 1) lies on a boundary. 

Recursion (|1.3() defines a non- autonomous dynamical system acting on the space 
of errors with a time-dependant parameter that belongs to the polytope. We are 
going to focus our attention on an associated dynamics in the ambient space of the 
polytope. We make the following changes of variables: 

(1.4) x(A:) = 7(/c) + e(fc - 1) , 
and 

(1.5) V{k) = v{x{k)) = v{j{k) + e{k - 1)) . 

Then adding j{k + 2) to both sides of 1(1. 3|l and reducing indices, we get 

(1.6) x{k + 1) = x{k) + 7(fc + 1) - v{x{k)) . 

Recall that v{x{k)) is chosen as the closest vertex to x{k) with some tie breaking 
rule. The orbit x{k) can be expressed in terms of the mapping t/)^ in 1(1.1(1 by 

(1.7) x{k) = x{k - 1) + (7(fc) - vix{k - 1))) = (f^^kMk - 1)) , 
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and bounding it in Euclidean norm is answered by Theorem ll.il Bounding e(fc) 
is equivalent to bounding x{k), and an immediate corollary of Theorem II. II is the 
following convergence of averages: 



1.1. The greedy algorithm in digital printing. To see the relevance of the 
above problem of approximating the 7(fc)'s by the V{k)^s to color printing requires 
some background in color theory. The color printing alluded to concerns digital 
halftoning which is the printing technique for imitating full color natural images 
(as opposed to colors placed at random on a page) by a checkerboard of dots limited 
to a few available colors. The simplest halftoning problem is that of obtaining gray- 
scale images with only black and white dots. 

Color theory itself is a fascinating subject which we can only briefly touch upon. 
It has long been known that our perception of color has a vector space model. 
Electromagnetic spectra are so rich that to describe them requires an infinite di- 
mensional function space. Nevertheless, experiments indicate that the human visual 
system reduces the visible part of this space to a three dimensional convex cone. In 
1931 the Commission Internationale de L'Eclairage (CIE) devised a certain three 
dimensional coordinate system called the tristimulus space which quantifies this 
cone. A complete account of how this coordinate system is contrived from color 
matching experiments is to be found in |30| . Our perception of light color is gov- 
erned by a set of rules known as Grassman's laws, the consequence of which is that 
two different linear combinations of color vectors appear as the same color perhaps 
with different radiances, if their sums are collinear. Consider a section of paper 
on which a color image is to be printed using a digital printer, typically a laser or 
ink jet printer which respectively uses toner or ink to render colors. We shall side 
step the nonmathematical practical difficulties involved in halftoning and assume 
an idealized situation: namely, the paper upon which an image is to be printed 
is partitioned into a checkerboard of tiny squares called pixels after a term used 
in color TV. We shall label pixel locations as we would the entries of a matrix. 
On each pixel is deposited a uniform color restricted to a very limited number of 
available choices which we shall call pixel colors. The problem is how should one 
choose pixel colors to render an acceptable full color image. 

Over a small area that still consists of a multitude of pixels the eye averages 
the received light. This means that we do not perceive colors on a page by taking 
positive linear combinations of color vectors, but rather convex combinations. The 
pixel colors commonly available to a printer are the standard inks or toner col- 
ors: cyan(C = (21,27,72)), magenta(Af = (33,18,22)), yellow(y = (65,76,14)). 
In addition red(i? = (30,18,7)), green(G = (11,22,13)), and blue(B = (9,7,20)) 
can be gotten by some kind of pairwise mixing of the above colors. Also avail- 
able are the ink or toner for black(iC = (5,6,6)), and white(VF = (84,87,105)) 
which is taken for the color of paper in the absence of any ink or toner. In paren- 
theses are typical tristimulus space coordinates for color vectors. The finite set 
{C, M, y, R, G, B, W, B} of ideally available printer color vectors are vertices 
of a convex polytope, in fact topologically a cube in three space. For a general 
reference about how color theory affects digital color printing see [O] . 



(1.8) 
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Since color usually varies with low spacial frequency over many portions of nat- 
ural images, it seems natural to simplify the halftoning problem by asking how 
should one choose pixel colors to represent a single uniform non-pixel color. Fur- 
thermore, since averaging plays such an essential role in color perception, it seems 
appropriate to recall that the uniform distribution theory developed by T. van 
Aardenne-Ehrenfest, H. Davenport, K.F. Roth and others. Their results concern- 
ing the discrepancy between average values and expected ones can be applied to 
the problem of approximating arbitrary uniform greys over two dimensional regions 
by combinations of black and white squares on a lattice. From a theorem of Roth 
(see 0, p. 3) one can deduce that it is not possible to bound the cumulative er- 
ror simultaneously for all grey levels and all rectangles. This is in contrast with 
printing in one dimension-i.e. printing on a line instead of a page- where it is 
well known, and as Equation Hl.l()|l below shows, that the cumulative error can 
be very well bounded on all segments, not only for grey levels but also for color. 
The discussion so far seems to suggest that, while there is no inherent difhculty 
in performing digital printing on a line, there may be one a page, where it is of 
practical importance. However, the issue is not one of bounding cumulative errors, 
but rather of limiting their growth rate. Furthermore, a two dimensional region R 
of pixels over which the eye averages is special. 

For the sake of simplification, let R he a. n x n square of pixels. The local 
cumulative error E{R) made over the region R is given by 

(1-9) E{R)^ (7(*,J)-^^(^,J))- 

In the terminology of printing, the 7(i,j)'s are called inputs and the y(i,j)'s are 
called outputs. The object of halftoning is to chose outputs so that the average 
cumulative error ||£'(i?)||/n^ is small. With respect to a polytope P in color space 
whose vertices are the inputs, as we have seen, the greedy algorithm of Equation 
(11.31) will achieve this. This is one of the methods of halftoning which is called error 
diffusion |10l 1111 12 71 1131 135| . We call this particular version simple error diffusion 
Since this involves a sequential process, pixels must linearly ordered, the usual 
order being the lexicographic one: namely, for an M x array the (j,j)*'* pixel 
gets index k — i ■ N + j where the indexing starts at 0. If _B is a bound on e{k), 
then 2B is a bound on the size of cumulative error 

L+n-l 

(1.10) e{L + n-l)-e{k)= ^ {lU)^V{j)) 

due a run of n pixels starting at pixel L.. Since there are n such runs in R, we have 
that ||i?(i?)|| < 2Bn which leads to small average size of cumulative error for large 
enough n. 

For two dimensional digital printing simple error diffusion is deficient because 
according to Equation 1)1. 3|l . even forgetting about the difficulty associated with 
changing lines, only pixel k — 1 has a direct influence on pixel k. Other nearby 
pixels have no direct influence. Even using more exotic orderings cannot eliminate 
this flaw. However, one can simultaneously improve the influence nearby pixels and 
deal with edge effects by using a general error diffusion which involves the greedy 
algorithm with respect to a weighting of past errors as follows. Let , 1 < i < m, 
be a system of weighting factors chosen in most cases for each fc to be a probability 
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Table 1. Table of nonzero weights 



vector. Indeed we shall always assume that the weighting factors are chosen this 
way: i.e., 



yf > 0, and = 1. 



Define 



(1.11) 



e{k + l) 



.1=1 



where V{k + 1) minimizes ||£(fc + 1)|| with some tie-breaking rule. In the lexi- 
cographic ordering of pixels, even though those to the right and below have no 
influence on the current pixel, this method of weighting the past orbit seems to be 
a much better approximation to the eye's method of averaging than that of giving 
all the weight to the pixel immediately preceding the current one. A typical scheme 
for weighting factors for pixels not near borders is shown in Table (^) where pi is a 
twelve-dimensional probability vector. Since this vector is allowed to depend on fc, 
adjustments can be made to deal with pixels near borders. 

The first such weighting scheme was introduced in "TI7 by R. Floyd and L. 
Steinberg in the seminal paper where they first described diffusion as a method for 
digital printing in 1975; this original scheme only involved pi, ps, and pe, all other 
weighting factors being 0. Two more elaborate schemes can be found in JI] and 
[ST] using the full Table and a lot more have been devised and utilized since. In 
these two schemes, as is often the case now for square pixels, the highest weights are 
given to pi = ps, the next higher to pi = pe and so on , with value decreasing with 
the distance to the fc"^ pixel being treated (for specific values and other information 
concerning digital halftoning see the above references as well as |T3 and Bll ) . 

For general error diffusion the relation oi e{k) to cumulative error is more com- 
plicated than Equation (|1.1U|) . Nevertheless we shall show that, as in the case of 
simple error diffusion, for general error diffusion ||_B(i?)||/n'^ 0. 

First we have the following corollary to Theorem ll.il 

Corollary 1.2. For general error diffusion ||e(fc)|| < B for some positive number 
B. 



Proof: Setting 
(1.12) x{k + \) 



Y,w'le{k + l-i) 



+ 7(^ + 1) 



(in the literature this is called the modified input)we have 



e{k) — x{k) — v{x{k)), 
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where v(x(k))=V(k) is the closest vertex to x(k) with some tie-breaking rule. Then 
(1.13) x{k + 1) = E [^(^ + 1 - - v{x{k))] + 7(fc + 1) 



1=1 



(1.14) 



The fact that x{k) of (|1.13|) is bounded, and equivalently that e(fc) is bounded, is 
a consequence of the convexity of Q in Theorem (|l.l|l . □ 

It is folklore knowledge in digital printing that for an n x n square R, the average 
error ||i?(i?)||/rt^ goes to zero with n (see for example [HIE])- For the sake of 
completeness we indicate a proof: more precisely that ||£'(i?)|| — 0{n). Suppose 
that the pixel locations of the nxn square R contained in an M x iV array, N >> n, 
are L + kK + j, < j,k < n — 1. Then from Equation (|1.11|) we have 

(n-l) L + feJV+n-1 

(1.15) E{R)= E E (7W-ni)) 

/c=0 j=L+kN 

m n-1 L+kN+n-1 

= E^* E E - eO' - 0) 

i=l k=0 j=L+kN 

The only nonzero weights are {wi, W2, WpN+q '■ 1 < P < 2, —2 < q < 2}. 

Changing the limits on the index j, and interchanging order of summation (|1.15|) 
we get 



n-l 



EiR) = J2^, E 



j=0 Lfe=o 



■ji-i 



n-l 



E e(j +L + kN)-J2 <J +L + kN-i)) 



k=0 



Let 
(1.16) 



E, 



E 

J=0 



n-l 



E e{j + L + kN)-Y^ e{j + L + kN-i) 



.k=Q 



k=0 



Then for < p < 2, -2 < g < 2, 

n—1 fn- 1 n—1 

(1.17) EpN+, = E E ^(J' + ^ + - E <^ +L+{k-p)N-q) 
j=a Lfc=o fc=o 

n—1 n—1 n— q— 1 n— p— 1 

= E E ^(^■+^+^^)- E E <j + L + kN). 

j—O k—Q — Q k——p 

We leave to the reader the subsequent calculations, the result of which is 

\\E{R)\\ = 0{n). 

1.2. A pursuit problem. The greedy algorithm also leads to a solution to the 
following pursuit problem. Suppose a predator is chasing a prey in a polytope, and 
the following conditions are imposed. The positions of the predator and prey at 
time t„ = i„_i + i are denoted by p{n) and q{n) respectively. We take n to be 
greater than zero and assume for that p(0) and g(0) are at two distinct corners, but 
any p(0) ^ q{0) would do. Between times tn and time tn+i the prey moves in the 
direction of any point -f{n+l) in the polytope with speed while the predator 
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moves at unit speed but restricted towards some vertex, say V{n + 1), the choice of 
which being dictated by a pursuit strategy. Notice that the prey has the advantage 
over the predator of freer movement but the disadvantage of slower speed (although 
less and less so). The rules of movement are expressed by the following formulas: 

(1.18) p{n + 1) = p{n) + -{V{n + 1) - p{n)) , 

n 

(1.19) q{n + 1) = q{n) + -^{j{n + 1) - q{n)) . 

71+1 

Defining a difference in positions per unit of time as 

e{n) — n{q{n) - p{n + 1)) , 
we have using 1)1.19(1 at time <„+2, and 11.18|l at time tn+i that: 

(1.20) e{7i + 1) = e{n) + j{n + 1) - V{7i + 2) . 

The reader might well anticipate that the pursuit strategy is the greedy algorithm: 
namely, the predator chooses vertex V{n+ 1) which minimizes of e{7i) with respect 
to the Euclidean norm, with some tie breaking rule. In terms of the notation 
previously introduced 

V{7i + 1) = v{e{7i ~1)+ 7(n)) = v{{n. - l){q{n - 1) - p(n)) + 7(71)) . 

Thus 

(1.21) e{n + 1) = e{7i) + j{n + 1) - v{e{7i) + 7(71 + 1)) , 

which we recognize as Equation II. 31 It then follows from Theorem 1 1.1 1 that 

\\q{n) -p{n +1)\\^0, 

which implies 

\\q{n)-pi7i)\\^0. 
In this sense the predator catches the prey. 



1.3. Instances of the greedy algorithm in classical mathematics. The greedy 
algorithm with which we are concerned in this paper leads one into a realm of 
tremendous richness. Consider the simplest input case where j{k) = 7 is constant. 
The performance of the greedy algorithm reduces to the study of iterations of a sin- 
gle map (j)-y. This map is a piecewise isometry. In general the analysis of piecewise 
isometrics can be of incredible depth and difficulty (see for example ^). Indeed, 
the simplest of all cases, the constant input case where the polytope P = [0, 1] the 
unit interval has significant mathematical substance. One easily checks that the 
interval / = [7—^,7+5] is invariant and absorbing under Both the interval 
/ and [0, 1] are fundamental regions for the action of Z on R and 

(f)-f{x) — X + J — v{x) ~ X + J mod Z. 

In other words, both / and [0, 1] can be identified with the circle R/Z and (f)^ with 
a rotation by an angle 7. This is one of the standard examples studied in ergodic 
theory and has the property of unique ergodicity when 7 is irrational. In Figure 1 
we have drawn the case in which < 7 < 1/2. A consequence of unique ergodicity 
is that for irrational 7 and any G /, the ergodic averages of the characteristic 



8 R. L. ADLER, B. KITCHENS, R. MARTENS, C. PUGH, M. SHUB, AND C. TRESSER 
Y-1/2 

1 

T 1/2 



Figure 1 . Diagram for simplest case 

fmiction X[i.-y+^]{^o + ^7 mod 1) converges to the measure of [^,7 + i] which of 
course is equal to 7. For this particular interval, H1.8|l and the equation 

X[i,^+i]ixo + mod 1) = X[i,7+i]('/'7(a;o) = vix{k)), 

where x(0) — xq, leads to the same result even for rational 7. 

The type of sequence v{x{k)) of O's and I's generated by the recursion 
is known as a Sturmian sequence. Such sequences were studied by Morse and 
Hedlund [2] ■ They have the property that there is a limiting frequency of 1 's and 
the distribution of O's and I's is as even as possible. Arbitrary Sturmian sequences 
can be identified with orbits of homeomorphisms of the circle, but only those that 
appear for rigid rotations get generated by the greedy algorithm in the constant 
input case (see |23 and (2H|). These sequences provide a symbolic counterpart to the 
theory of Diophantine approximation (see |2H! and ^S]) and have a long history. 
Johan Bernoulli was apparently the first to discover what came to be known as 
Sturmian sequences. He did that to describe an easier method of interpolation in 
astronomical tables Q; Christoffel applied such sequences to modular arithmetic 
121 ; Smith used them to study well distributed sequences of two symbols ■ Both 
Christoffel and Smith seemed to have been unaware of the earlier work of Bernoulli. 
Bernoulli's work was later revisited by Markov ^H]- A generalization of Smith's 
view of the problem was considered in where the circle was replaced by arbitrary 
dimension tori, and the corresponding multidimensional continued fractions were 
described. 

One more feature of the constant input case is the fact that under the action of 
(j)^ on R the interval I is invariant and absorbing-i.e., every orbit eventually enters 
and stays in I. We expect that the following generalization holds in all dimensions, 
the proof of which we shall defer to a later work. For any lattice C with a preferred 
basis in R^, let P be the standard fundamental parallelepiped in the basis. Then 
for any 7 € P there exists another fundamental region Qj which is an invariant 
absorbing set under the action of (j)-^; and on Q-y the map 0^ is a rotation by 7 on 
R^/£. The same results holds if P is a simplex instead of a parallelepiped leading 
however to a different fundamental region Q^. Furthermore, there is a convergence 
of ergodic averages for special sets in the simplex case similar to the interval case. 
Finally, it is conceivable that this result could lead to higher dimensional analogues 
to Sturmian sequences. 

1.4. Outline of the paper. In §2 we describe how the proof of Theorem 11.11 
evolved from dimension 2 to the general case. In §3 we prove Theorem 11.11 based 
on two propositions. Proposition 13.11 and Proposition 13.21 which we prove in §4. 
In §4, we show that the proof of Proposition 13.21 itself easily follows from a third 
proposition, Proposition l4.1l Proposition l4. ll posits the existence of a special subset 
of the {N — l)-dimensional sphere upon which the construction of the invariant set 
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Q is based. We postponed the proof of that Proposition to §5. Item (3) of that 
proposition is where aU the difhculties and subtleties of our proof of Theorem 11.11 
he. FinaUy, at the end of §5 we have finished assembling all the elements for a 
complete proof of Theorem 11.11 and in §6 we present two Theorems concerning 
some additional properties which are byproducts of our method of proof. 



2. How THE PROOF EVOLVED FROM DIMENSION 2 TO THE GENERAL CASE 

We shall adhere to the convention of writing the inner product of vectors as a 
dot product, the Euclidean norm as || • ||, and Euclidean distance between a pair of 
points or between a point and a set as •) . 



Definition 2.1. Given a set of points wq, . . . , vm-i in any dimension, the Vorono'i 
region Ry. of Vi is defined as 



Rvi = f^{x ■■ \ \x ~ Vi\\ < \ \x - 



hence a polyhedron being the finite intersection of half-spaces. It is the set of 
points closer to Vi than to any other point in the set (or possibly no further than 
equidistant to some). Also let v{x) denote the closest vertex to x with a tie breaking 
rule such as choosing the vertex of smallest index in case of ties. 

Now let P be a convex polygon with vertices vq, . . . , vm-i ■ To make the definition 
of the map (f>-y{x) more specific, the tie breaking rule for v{x) will be for us the 
smallest index i such that x G R^-. For example, if P is an interval [woj^i], then 
the Voronoi regions are: 

Ry^ = (-00, J , 



and 



Rv^ = ( — — 



The mapping : R ^ R is given by: 

4>-y{x) = a; + 7 — v{x), where v(x) 



Theorem 12 . II and Theorem 12 . 21 are respectively dimension 1 and 2 (intervals and 
polygons) versions of Theorem (|l.l|l . Proofs of them appeared in 

Theorem 2.1. Let Qt denote the interval generated by moving the ends of P — 
[vq-, Vi] outward a distance t > : i.e., Qt = [vq — t^vi + t]. Then for any t > 
and 7 e P 

(f>j{Qt) c Qt. 



Let Vi, i = 0, . . . , M— 1, n > 2 indexed in clockwise order P and Uj, j = 1, . . . ,n 
the unit normal vectors to edges Vj-iVj of P. Then Uj ■ {vj — Wj-i) = 0, and we 
can write rij ■ vj-i = Uj ■ vj = dj. 



10 R. L. ADLER, B. KITCHENS, R. MARTENS, C. PUGH, M. SHUB, AND C. TRESSER 





Figure 2. Polygons P, Pt, and Voronoi region R 



Theorem 2.2. Let Qt denote the polygon generated by moving the edges of P 
f]-{x:x-nj<dj} perpendicularly outward a distance t > : i.e., 



Qt = f^{x : X ■ Hj < dj + t}. 



See Figure\^ There exists T such that for any t > T and 7 G P 

(pfiQt) c Qt. 



The above method of constructing an invariant Q by outwardly translating the 
planes of faces of a polytope P, even by different distances, is doomed to failure in 
dimension A'^ = 3. as shown by the following counterexample. 

Consider the octahedral polytope abcdef ghijkl in Figure The faces abcdef 
and ghijkl are parallel to the plane (0, 0, 0)(0, 1, 1)(1, 0, 1) cutting inside the cube as 
shown and equidistant from the points (0, 0, 0) and (1, 1, 1) respectively. The face of 
the Voronoi regions Rg separating g from h is in the plane perpendicular to the top 
face of the cube through the midpoints of segments do and kg. Invariance fails at the 
point the midpoint m of do because here the vector d — g sticks out of the polytope. 
No matter how far the top plane of the cube is translated upward this vector will 
still stick out the resulting polytope at such a point. So this means that the faces 
abcdef and ghijkl must be translated far enough apart. But then a new difficulty 
appears: the face of Voronoi' region Rc lies in the plane determined by (0,1,1), 
(1, 0, 1), and the midpoint of be. The vertex (1, 0, 1) of the cube lies in this face, but 
here the vector b — c sticks out of this new candidate for an invariant polytope. On 
one hand this second difficulty cannot be overcome by moving the front and side 
faces outward, but on the other hand moving the top face up reintroduces the first 
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Figure 3. Polytope Counterexample 



difRculty. We leave it to the reader to verify that one problem cannot be overcome 
without introducing the other: for details and the general case of dimension N > 2, 
see 

So another idea for a general construction of an invariant set is needed (special 
constructions for special polytopes will be discussed in JHI)- Imagine the limiting 
case of a sphere centered at the origin scaled down to the unit sphere as the radius 
goes to infinity. The edges at a vertex are normals to hyperplanes through the 
origin, which determine the limiting Voronoi' region at that vertex. The following, 
which we shall verify later for all dimensions, is easy to see in the two dimensional 
case below. It is the most important element of the proof. The cone of normals to 
the supporting hyperplanes at a vertex agrees with the Voronoi' region at infinity of 
that vertex (See proof of Proposition l3.1|l . Thus the tangent space of the sphere at a 
normal in the Voronoi region of a vertex is a supporting hyperplane of the polytope 
at that vertex. In this idealized situation the polytope has shrunk to a point at 
the origin; and the unit sphere (equivalently the sphere at infinity) is invariant. 
Backing off from infinity to a sufficiently large sphere, the invariance fails only on 
neighborhoods of Voronoi boundaries. Resolving this difRculty then becomes the 
main task of the proof. To do this we construct a special subset of the unit sphere 
the tangent planes of which are supporting hyperplanes of a proto-invariant convex 
set Qao- Then the set Q = pQoo for large enough p will be the desired invariant 
convex set. 

Before launching into the formal proof for all dimensions, we illustrate the idea 
for some low dimensional polytopes indicating the induction upon which the general 
construction is based. 

First we take up a construction of Of Q viala setil in dimension 2. Consider 
the idealized limiting circle scaled down to the unit circle as depicted in Figure 
Q). Let vi,vq,V2 be three successive vertices of a polygon, and vi = vi — vq and 

= i'2 — wq. The points b and g lie the intersection of rays perpendicular to vi 
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Figure 4. Construction of Q in diniension2 



and 1^2 with S. These rays are parallel to the unbounded edges of the of the Voronoi 
region Ry^ for the vertex vq. Using the same letters to denote corresponding points 
on a dilated circle pS, the intersection of the infinite bounding edges of the actual 
Voronoi region with the circumference of the dilated circle for large p would be in 
the arc ad and in the arc eh. The only place a vector 7 — wo would stick out of the 
circle for some 7 G P would be at points b and g {e.g., j — Vi at 6) and nearby 
points relatively speaking. So in order to avoid this difficulty, we remove small open 
arcs ad and eh of S, retaining the points b and g, and carry out this procedure for 
each Voronoi region Ry. . The size of the arcs is subject to the restriction that the 
removed arcs don't interfere with one another. It can then be seen that, if Qoo is 
the convex set supported by tangent lines to the unit sphere at points of CI, then 
the convex set pQoc is invariant for large enough p. 

A construction for Q via a set fiin dimension 3 is similar. For the sake of 
simplicity, suppose that three edges meet at vertex vq and that vi, V2, and V3 
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Let vi — vi — vq, 





Figure 5. Construction of Q in dimension 3 

1^2 — V2 — V0, and I's — V3 — vq be the normals to the faces of the the Voronoi' regions 
Ryg . The planes defined by these normals intersect the idealized scaled limiting unit 
sphere in the spherical triangle abc. To construct f2 we remove from the unit sphere 
small open spherical caps centered at a, b, and c, but retain these centers. In the 
complement of what has been removed, we remove open strips centered on the 
edges of the spherical triangle, but retain their center lines which in this case are 
the closed arcs c?e, fg, and hk. In Figure (0) we have illustrated the construction 
for Voronoi region . We carry out this procedure for each Voronoi region. Again 
caps and strips should be chosen so that they do not interfere with one another. 
The set Vt is what survives. As before, it can be seen that the failure of invariance 
of large spheres is overcome by the set ^2, and the convex set Qoo supported by 
tangent planes to the unit sphere at points of fi, whereupon Q = rQoo is invariant 
for large enough r. 
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Note that the boundary of the removed caps from the 2-sphere in 3 dimensions 
has the same structure as in 2 dimensions. Also the boundary of a cross section 
of a removed strip from the sphere has the same structure as a hne through a 
boundary of a removed arc on the circle (like the situation on the arc ad). Thus 
we get a glimpse of the induction about to be carried out. 

In dimension 4, the intersection of a Voronoi' region and the 3-sphere is a spherical 
polytope. To construct f2 we proceed exactly as before. We remove small open 4- 
spheres about the vertices of the polytope retaining the centers. From what remains, 
we remove open tubes around the edges, retaining the edges. Finally from what 
still remains of the unit sphere we remove slabs about the faces of the polytope 
while retaining the faces. We do this for each Voronoi region. Again a restriction 
on the removed sets is that they don't interfere with one another. 

Again we note that the boundary of a removed spherical neighborhood from the 
3-sphere in 4 dimension is a 2-sphere having the same structure as the 3-dimensional 
n. The retained center of such a removed spherical neighborhood is not visible in 
the 3-dimensional subspace in which the bounding 2-sphere lives. Rather the center 
of this 2-sphere is a projection of the retained point in f2. This corresponds to what 
is happening in the previous dimension. In Figure ((SJ the point a is not in the plane 
which contains the boundary circle of the removed spherical cap on the 2-sphere. 
The center of this circle is the projection of the point a to this plane. Going down 
one more dimension, this corresponds to the fact that b does not lie on the line 
through a and d. A cross section of a removed tube has the same structure as the 
boundary of a removed cap in dimensional 3 which in turn has the same structure 
as the removed circle in dimension 2. 

Remark. Returning to Figure we note more points could have been removed 
between d and e without compromising invariance. Similarly, the same can be said 
for removing more points from Q in any higher dimension. Our actual construction 
of Q will take advantage of this fact, and we shall remove much more from the unit 
sphere than what we have just described. This has two virtues. First describing 
removed sets is slightly easier. Second it leads to a simplified proof from which we 
get for free Theorem 16. 21 

For the construction in dimension 2 of this more truncated fl, let p and p' be the 
points of intersection of the unit sphere with the line normal to v = Vi — Vj for a 
pair of vertices Vi and Vj . Remove from the unit sphere small open arcs (symmetric) 
about p and p' but retaining p and p'. Do this for every pair of vertices v^ and Vj 
of the polytope. The only restriction is that removed arcs be small enough not to 
interfere with each other. The set f2 is smaller than in the previous construction 
which results in fewer supporting lines for Qo^ which then results in a larger Qoo- 
Nevertheless, the ultimate Q will still be invariant. 

To construct fl in dimension 3, consider all possible intersections with the unit 
sphere of planes whose normals are given by v^ — yj for every distinct pair of 
adjacent vertices v^ ^ Vj of the polytope P. These intersections consist of one 
dimensional lines, and two dimensional planes. We first remove caps centered at all 
the points that are intersections of the one dimensional lines with the unit sphere 
while retaining the center points. From what is left we remove all strips centered 
at great circle arcs that are intersections of two dimensional planes with the unit 
sphere while retaining the central arcs. Once again the removed sets must be at 
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least small enough so that they do not interfere with one another. To facilitate our 
proof it turns out to be convenient to make these removed sets even smaller than 
just the requirements for non interference. As before, the resulting Q leads to an 
invariant Q. 

The above will be the procedure for constructing Q in all dimensions. Showing 
that a supporting hyperplane tangent to the unit sphere at a point in Q is not cut 
off at the wrong place by another one is delicate. This is clear in dimension 2, but 
far from obvious in higher dimensions. 

3. Proof of Theorem 

We shall use freely some notations, definitions, and basic results from classical 
convex geometry (see for instance |26) or |27j and references therein). 

3.1. Notation. 

• P is a polytope-namely, the convex hull of a finite set (vertices) {vi, . . . ,vm} 
of extreme points (vertices)in M^. By the classical Finite basis theorem for 
polytopes attributed to Minkowski [1896], Steinitz [1916], and Weyl [1935] 
(see 123; P- 89), P is a bounded polyhedron. 

• C{vi) denotes the cone of outward normals to the supporting hyperplanes of 
P at vertex u^; and Ci{vi) the unit normals in C{vi). See spherical triangle 
abc in Figure 10). 

• As a polyhedron, P can be expressed by 



3.2. Proof of the main result. We shall assume two propositions, the proofs of 
which shall be deferred to the next two sections. 

Proposition 3.1. For each e > there exists D > such that ?/ ||a:;|| > D then 



Proposition 3.2. There exists a closed subset fl of the unit sphere in R centered 
at the origin with the following properties: 

(1) n is not confined to any half-sphere; 

(2) there exists ane, < e < 1, such that if y € Ci{v) and lu G fi — Ci(u) then 
there exists uJ £ Ci{v) n fl such that 



(3.1) 




all Vi ujeC{vi) 



(3.2) 



d{uj',y) < d{uj,y) ~e. 



Proof of Theorem 11.11 



Let 



(3.3) 




wen 
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and 

(3.4) pQoo^ f]{x\x-LU < p}. 

wen 

We shall show that for all sufficiently large p the set pQ oo satisfies Theorem ll.il It 
is clear by definition that Qoo and hence pQoo are convex. If Qoo were not bounded, 
it would contain a ray R+x = {tx\ t > 0}, and the half-sphere determined by the 
hyper-plane normal to x would contain fJ, contradicting Proposition l3.2l Thus Qoc 
and pQx> are bounded. 

We must show that there exists an R such that 

(3.5) u) ■ 4>-y{x) = w • (x -f 7 — v{x) < p 
for all X e pQoo, p > R, 7 G P, and w G fi. 

Case I: a; e C{v{x)). 

By (|3.1|l and (|3.4|l . for 7 e P we have 

(3.6) ll) ■ X + Lu ■ ~ v{x)) < p . 

The inequality is satisfied no matter what the size of ||x|| and p > are. 

Case II: to ^ C{v{x)). 

Let e be the one given in Proposition 13. 21 From Proposition 13. II there exists Rq 
such that for ||a;|| > Rq there exists y £ Ci{v{x)) satisfying 

rf(^,y)) <e/4, 

By Proposition 13 . 21 there exists w' G Ci{v{x)) n such that 

d{w,y) > d{uj',y) + e. 

By the triangle inequality 

d(..,^)>d(a;,y)-%,^) 

> diuj',y)+s-d{y,j^^) 
>rf(^,c.')+e/2. 

Let 6 be the angle between x and uj, and 9' between x and lo'. Using the fact that 
differences in arc lengths are greater than differences in corresponding chord lengths 
(an exercise-hint: make chords parallel) we get 

(3.7) e_0'>d(^,^)_d(^,^')>£/2 
for the arc length difference. Thus 

(3.8) cos(6') < cos(6'' + e/2) < cos{e') cos(e/2) , 
and we arrive at 

X f X , > 

UJ ■ < uj ■ T^fl - e) 
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where e = 1 — cos |. Note < e < 1, depending only on e. We then get 

(3.9) uj-x + LJ- {j- v{x)) < J ■ x{\ - e) + diam(P) . 

We can therefore choose Ri > Rq such that for ||a;|| > Ri the right hand side of 
(Ol) is < p. Let 

(3.10) R = Ri + diam(P) , 

and let x G pQoo where p > R. By what we have just shown, if ||x|| > Ri then 
X satisfies equation (13.51) : so ^^(x) G pQoo- Finally, if ||a;|| < Ri then x G pQoo ■ 
because \\x — (j)^{x)\\ < diam(P). □ 

Remark 3.1. At this point we have also proved that given any compact set K the 
set Q can be chosen to contain to contain K. 

4. Proofs of Propositions H3.1|l and (|3.2|) 

4.1. Proof of Proposition 13. fl By a standard argument in linear programming- 
namely, for each vector u the function u ■ x achieves a maximum over P at some 
vertex Vi which means that u-{x — Vi) < for x G P, or equivalently that u S C'{vi)- 
we have = [J,C{vi). 

Furthermore, by the Decomposition Theorem of Polyhedra due to Motzkin [1936] 
(see (27], p. 88) a Voronoi region P„;, being a polyhedron, can be written 

where P^. is a polytope (recall a polytope is a bounded set) and K^. is a polyhedral 
cone. Since = IJ - Ry^ , it is easy to prove that R^ = IJ ■ Ky. as well. 

We show that K^^ = C{vi) as follows. Let u G C{vi), and t > 0. The vector u 
determines a supporting hyperplane for P at Vi which implies that Vi is the closest 
vector of P to the vector tu + Vi : so tu + Vi £ Ry. . Thus from the polyhedral 
decomposition theorem tu + Vi — p{t) + k{t) for all i > where p{t) e Py^ and 
k{t) e Ky^. Therefore, 

(4.1) ^^p{t)-v. 
^ ^ t t 

The first term on the left in equation 14.1|l converges to as i — > oo so that the 
second term, which always belongs to K{vi) by convexity, is bounded for all i > 0. 
Thus there is a convergent subsequence whose limit is in K{vi) which implies that 
u e Ky. : hence C{vi) C Ky^. 

Since the two collections of C(ui)'s and of Ky.''s form partitions of M^, the 
reverse inclusion follows easily from the fact that the interiors of C{vi) and C{vj) 
are disjoint for i ^ j, and similarly for the interiors of Ky. and Ky^ . 

Given any x G R^, we have just shown 

X e Py(x) + C{v(x)) , 

so there exists p{x) £ Py(x) such that x — p{x) G C{v{x)). Theorem 13.11 now follows 
from the easily proved limit 

.. X X — p(x) 



\x\\ ||a;-p(a;)|| 



, as ||a;|| — > oo . □ 
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4.2. Proof of Proposition 13721 The proof of Proposition 21 follows easily from 
the following proposition, the proof of which has been postponed until the last 
section. As previously stated, this proposition is the heart of the matter 

Proposition 4.1. There exist a partition of each Ci{vi) into a finite number of 
connected sets a and a closed subset il of the unit sphere S in M.^ centered at the 
origin with the following properties: 

(1) n is not confined to any half-sphere; 

(2) for each a C C\(vi) there is a eg- > such that 

where \J^{A) = {x\ d{x,A) < e} denotes the e— neighborhood of a set A; 

(3) the closest points uJ€^toy'Ed'(l Ci{vi) are in the same set closure a. 
In particular ct H 17 7^ 0. 

The sets a € Ci(vi) will be chosen as cells of a simplicial decomposition induced 
by great spheres on S as specified in §5. 

Proof of Proposition l372l 

We must prove that there exists e such that for all Ui and all y e Ci{vi) 

(4.2) d{y, n \ Ci{v,)) - d{y, Ci(f,) nn)>e. 

For cr C Ci{vi) define a function /g. : (i ^ M by 

fg{y) = d{y, n \ (a)) - Ci (v,) n n) 

where Eo- is given by Proposition 14.11 Then /o- is continuous and strictly positive 
on a compact set: hence assumes a minimum e^ > 0- We claim e — niineo-, where 
the minimum is taken over all a C Ci{vi) and all w^'s, satisfies (|4.2|l . Observe that 

(4.3) diy, n \ Ci {v,)) - d{y, C\ (v,) nn)> d{y, d{y, Ci [vi) n , 

where y £ a C Ci{vi). Now by (2) and (3) of Proposition HTI the right hand side of 
63 is equal to d{y, n \ Ue„ (a)) - d{y, Ci{v,) n f^) = /^(y) > e. □ 

5. Proof of Proposition 14. II 

5.1. Notation. 

• S = S{c,r) C M"+^ ~ denotes an n-dimensional sphere with radius r 
and center c. We shall call the intersection of S with any fc-dimensional 
affine subspace of M"^^ containing c a great sphere of dimension fc — 1 or 
simply a great sphere. 

• For A a closed subset of S, the set Clos(a;, A) of points in A closest to a 
point a; G S" is given by 

Clos(2;, A) = {a e A\yb e A, \a-x\<\b- x\} . 

• V = {Vi,V2, ■ ■ ■ ,Vs} denotes a collection of codimension 1 affine subspaces 
of M.^ containing the center of S. For any x in S, Vx will stand for the 
collection of Vi's that contain x. 
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For each a a subset of {1, 2, . . . , s} we define the subspace 

14 = n 

and the great sphere in S 

s-^ = 5 n v„ . 

By usual convention, Vg = and 5*0 = S. 

We call the pair S = (5, V) a stratification of S induced by the subsets of 
V. We call the lowest dimensional Sa containing a point x the stratum of 
X and denote it by Sa{x). It is well-defined, since the intersection of any 
such two strata would be a lower dimensional one. A stratification induces 
a partition of S into cells of dimensions 0, 1, ... , dim 5*, the cell containing 
X, denoted cr{x), being the connected component of {y g Sa{x)\ Sa{y) = 
Sa{x)} containing x. Two stratifications §i — (5i,Vi) and S2 = ('S'2,V2) 
are called isometric if there exists a bijection j : Vi — > V2 and an orthogonal 
afhne map O : Si S2 such that 

o{SinVi) = S2nj{Vi) 

for every 1^ G Vi. 

Let C be the afline subspace orthogonal to Va C passing 

through the center of S and V^{e) = Ue(c) n V^. 
The e-tubular neighborhood of Sa is defined as 

T,(5«) = (V«x V;^(e))n5. 

When the dim^Q, = k we also call T^{Sa) a k-tube. 

For every z e we denote the great sphere orthogonal to Sa at z by 
S^{z). Notice that 

- (R- (z-c) X V^)nS. 

For every z ^ Sa it is easy to check that the set 

Fa{z, e) = (M+ • (z - c) X V^{e)) D S . 

is a spherical cap. We call it an (a, e) — spherical cap with top at z or simply 
a spherical cap. 

For every z G Sa we define a neighborhood fibered by {a, e)— spherical caps 

Paiz,e)= y Fa{y,e), 

called the (a, e)-patch centered at z where a is given by the stratum con- 
taining z. 

The boundary dreiFa{z, e) of the spherical cap Fa{z, e) relative to the sub- 
space Sa{z) is given by 

drelFaiz, c) = (M+ ■ [z - c) X (5U,(c) H V^)) H S . 

The relative boundaries of spherical caps are Euclidean spheres though not 
great ones: rather they are analogous to meridians with the north pole at 
z. Note that 

dimdreiFa{z, e) = n — dim^a — 1 = dim5 — dim 14: • 
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• The stratification § induces the stratification {dreiFa{z, e), Vz) on the sphere 

drelFa{z,e). 

We leave to the reader to check the following relations: 

(5.4) T,{S^)= y F„{z,e), 

otherwise speaking, the {a, e)-spherical caps foliate the e-tubular neighborhoods of 

(5.5) F„(^,e)cU(^),(0); 

(5.6) Pc«(a;,e) C U(i_^^)^(a;); 

(5.7) S= \J S^{z); 
and 



dreiFc,{z,e) = S{c+ ^''' '\ z - c), e) n S^{z) . 
r 

We also leave as an exercise the following. 

Lemma 5.1. (Spherical Cap Charactization) A spherical cap is the intersec- 
tion of a half-space and a sphere, and conversely. 

5.2. Spherical caps, the Set n, and the Main Lemma. 

Lemma 5.2. (Spherical Cap Boundary Isometry) // = Vy and x,y G 

Sa{x), then {dreiFc,{x,e),Vx) and {dreiFo,{y,e),Vy) are isometric stratifications for 
any e > 0. 

Proof: Consider the plane oi x,y and c. This plane is a subspace of V^. Then 
the afline orthogonal map O, which rotates this plane about c rotating y to a; and 
which is the identity on the space orthogonal the plane, is the required isometry . 
□ 

Lemma 5.3. (Spherical Cap Intersection) Let Sa' C Sa- Given ^ < r, if 
X € Sa, and z € Sa', then either 

S^{x)C^Fa'{z,ii)=%, 

or 

(5.8) S^{x)^Fa>{z,ii) = Fa{x,8), 
where 5 and fi are related by 

(5.9) r^-i? = (r^ - 5^) (^ ^^"^ 
with r standing for the radius of the sphere S centered at c. 



Moreover, if ^ < \/2\/2- 2 ~ .911, and^ <t then 



(5.10) d{x,drelFa'{z,fJ,)) < tS . 
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Figure 6. Spherical cap intersecting a great 2-sphere 

Proof: Observe that if x ^ Fa{x,fj,) then S^{x) f] Fa'{z,^) ~ 0. Now let 
X e Fa{x,^j,) \dreiFa'{z,^) and X ^ z. then S^{x) r]Fa>{z,ii) ^ 0. Then S^{x) f] 
F„-(z,Ai) 0. By LemmaEl\S^(x]nF„Az. u) is a shperical cap and we are justified 
in denoted it by Fa{x, 6). We must now determine 6. 

Choose V £ dFa'{z,ii). Then v e dreiFa{x,6) n dreiFa'{z, n). Then the vectors 
v — c, x—c, z—c determine the three dimensional c-centered sphere S3 that contains 
the points v,x, and z. We assume that c is at the origin and the point z is at the 
north pole. The intersection of Fa' {z, fj,) with 6*3 is a spherical cap whose boundary 
is a meridian which is also dreiFa' (2:, /x) ft ^3; this meridian is determined by ji and 
the intersection of dreiFa{x,6) with S3. This intersection is an arc of a great circle 
whose endpoints lie on this meridian. Let y be the intersection of the great circle 
determined by x and z with this meridian. The projections of u on a; and z are 
P— (7^)2; and q = {^)z respectively. Note that 

(5.11) ^-./W^, 
and that 

(5.12) ^ = . 
We can write 

(5.13) v—( — :^)x + U, 

which exhibits the contribution u G C to v and yields 

(5.14) V ■ z ^ i—Y^)^ ■ ^ ■ 

Equation (|5.9|l then follows by substituting I5.11|l and (|5.12|) in H5.14|l . 

For the case x = z it is apparent that (|5.9() holds with 5 = ji. 

Finally, for the case x £ dreiFa' {z, fi) it is easy to see that (|5.9ll holds with S — 0. 

To obtain Equation H5.1()|l we first prove that y is the closest vector to x in 
dreiFa' {z, jj). Let V play the role of any vector in dreiFa' {z, /u), not just one of those 
at a distance 6 from the axis of S through x, and observe that d{x, v) > d{x, y). 
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From the fact that angle xpy is acute~i.e., 

(5.15) dix,y)^ <d{y,p)^ + dip,x)\ 

and the fact that the shortest distance to the circumference is along a radius-«.e., 

(5.16) d{p,x) <d{p,y), 
we get 

(5.17) dix,y)<V2d{y,p). 
Furthermore, from 

diy,p) = d{y,q) - d{p,q) , 
and the fact that Avpq is a right triangle and we have that 

(5.18) d(x, y) < V2{fi - - S^) 




A straight forward computation shows - < \/ 2\/2 — 2 is equivalent to 



(5.19) 
Thus 




For the case a; = z it is apparent that (|5.8|) holds with 5 = fi. 

Finally, for the case x £ dreiFa'{z, fi) it is easy to see that (|5.9|l holds with 5 — Q. 
□ 

The following lemma is an immediate corollary of the first part of the preceding 
one where z is replaced hy z' , x is replaced by is replaced by e', and e is taken 
to be a number not necessarily equal to 5. In fact Case I corresponds to e > i5. Case 
II to e = (5, and Case III to cither e < 5 oi an empty intersection. 

Lemma 5.4. Let Sa' C Sa- Given z e Sen o.i^d z' e Sa' ^ one of the following 
mutually exclusive relations hold: 

I. S^{z)nF^,iz',e')DF„{z,ey, 

II. S^{z)nF^,{z',e')=F^{z,ey, 

III. S^{z)nF^,iz',e')CF^{z,e). 

Remark 5.1. Observe that by the Pythagorean Theorem d{z,dreiFa'{z, fi))) > fi. 

The authors found schematic representations of Figure and Figure useful 
mental aids for arguments to follow. 
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Figure 7. Schematic aid 



Definition 5.1. (The set fl) With respect to an n-tuple e = (eo, ei, • • • , Cn-i) of 
positive numbers set 

dim Sa=j 

for j = 0, 1, . . . , n — 1. Next set 



f2o = ( U Sa) 



dim S„=0 



and 



Ofc = 0fc(eo,...,efc-i) = ( U Sa) \\jTjiej). 

dim Sa—fc j<k 

for fc = 1, . . . , n. Note that sets in this hierarchy may be empty because there are 
not enough afBne subspaces in V or the ej are too big. Finally, set 



k<n 



Definition 5.2. An n-tuple e = (eo, ei, . . . , CdimS-i) of positive numbers is said to 
decrease sufficiently rapidly and prevent improper tube interference with respect to 
a stratification {S, V) if 

(1) for A: > 0, ek/ek-i < t, with eo < rr, where r = -^^^ — , 

(2) for < fc < dim5-l andw e Qk, if C ^ V^,, then U8eJa;)nTe^(S'„) = 0. 

Item (1) will ensure that the e^'s decrease fast enough. The value of r is chosen 

so that 



(5.20) 



U8e.(^)CU 



{z). 



This relation implies that a 8e/;-neighborhood of any point in a sphere of radius 

Sk-i is contained in an octant of that sphere centered aroimd the given point. In 
subsection 5.4.2, claim 2, of the proof of the Main Lemma we shall need that the 
hypotenuse of a geodesic spherical right triangle is longer than the other sides. This 
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will be the case if the triangle lies in a neighborhood contained in an octant. In 
addition this value of r is smaller than the quantity \/ 2\f2 — 2 which appears in 
Lemma [5.31 It is also small enough so that each time a geodesic belongs to the 
discussion {e.g. see subsection 5.4.3, claim 3 in the proof of the Main Lemma) it is 
unique. 

Item (2), which might impose that ratios of successive are much smaller than 
r, guarantees that spherical caps fibering tubes in general do not interfere in unex- 
pected ways with spherical caps centered on points of Q and leaves enough room to 
simplify some arguments. We call attention to the reader that the subscript k in 
T^kiSa) is not necessarily the same dimension as that of Sa as it would be in the 
definition of Tfe in Definition 15. II 

Lemma 5.5. With respect to the stratification (S, V) let e = [sq, ei, . . . , e„_i) de- 
crease sufficiently rapidly and prevent improper tube interference. Then for each 
X G ljj<n there exist a largest integer < k{x) < n — 1, a great sphere Sa with 
dim^a — k{x), and a point z € Sa such that the {a, e^^)) -spherical cap Fa{z, e^^)) 
contains x and satisfies 

dreiFaiz, efe(,)) H (T_i U Ti U . . . Tfe(,)_i) = , 

where T_i = 0. 

Moreover, this {a, ek{x))~ spherical cap with top at z and containing x is unique. 
It will be denoted by Fa(x){z{x), ek(x)) to show all the dependencies on x, however, 
for the sake of brevity we shall also denote it simply by Fa{z, ek(x))- In the abbrevi- 
ated form it is to be understood that the dependence of a and z on x is adequately 
indicated by the subscript k{x) on e. Here the unique great sphere Sa containing z 
determined by x we denote by Sa(x)- 

Proof: Take the largest ki < n such that x E T^^ . By H5.4|l x belongs to some 
spherical cap Fa^izi, Ck^) where dimS'aj = fci. If ^rei-pQi (2:1, efc^ ) n Tj =0 for all 
j < fci we are done and the spherical cap Fa-^ {zi, efcj is unique. 

The argument for the uniqueness is the following. There are two cases: 

(1) zi e rifc, n Te,^ (5„ J ; 

(2) zi € Te^{Sa'),j < fci, dim 5a' ^ j, and Sa' C Sa^- 

In case (1) Fa-^ (zi, Cfci) C Use^^ (-^i)- The uniqueness follows from the fact that there 
are no improper tube interferences. 

In case (2) we have zi G T^j{Sa') with j < fci, say zi € Fa>{z,ej). Lemma [5.31 
implies that d{zi , w) < re^-^ for some ui G ftki (in that lemma the role of x is played 
by zi, S by efc^, ^ by e-j, and z by itself. Here FQi(zi,efcJ C Vsck^i^)- 

Hence also in case (2), the uniqueness follows by the way ek^ is chosen so that 
Faj(zi,efci) does not intersect any other (e^^ , a")-spherical cap with dim5Q" = fci. 

If dreiFai (zi , Cfci ) H T j 7^ for some fc2 < fci , take the largest such fc2 . By Lemma 
15.41 and definition of tki we get that x G -Fqi (zi, e^j-J C T^^ so that x belongs to 
a spherical cap Fa.^ (z2 , ) where dim Sa2 = ^2 ■ Either we are done and this cap 
is unique or we must repeat the procedure finding k^ and so forth. This procedure 
terminates at some fc > which defines k{x) and a unique spherical cap Fa{z, tk(x)) 
satisfying the conditions of this lemma. □ 



CONVEX DYNAMICS AND APPLICATIONS 



25 



Main Lemma 5.6. Let (5*, V) be a stratification with dim 5 = n and ri(e) = 
r2(eo, • ■ • , Cn-i) as in Definition \5.1\ If e = (ep, ei, . . . , e„_i) decreases sufficiently 
rapidly and prevents improper tube interference, then 



clos(x, fl) C a{x) . 

Before giving a proof of this lemma, we need some more notation. 

5.3. Induced Stratification and Related Constructs. Choose some a and set 
dim Va = fc + 1 , hence dim Sa — k . Pick any x Cz Sa, and set S — drei Fa ix,ek) and 
V = Vx ~ {Vi\ i G a} . Recall (see (|5.3f) ) that the sphere S has dimension n — k — 1. 
We denote the stratification induced by V on dreiFa{x,ek) by S = (5*, V). Let 

(5.21) e = (eo, . . . ,e„_A;-2) , 

For a C a we denote the great sphere 5 H Va by Sq. and the e-tubular neighborhood 
iV& X y/(e)) n of 5a in 5 by f,{S&) . Notice that 

(5.22) s& = s&ns, 

and that by H5.1|) 

(5.23) t{S&)^T,{S&)nS. 

Furthermore, one can check that the cells (j(x) of the induced stratification are 
obtained by intersecting the original cells with S: i.e., 

(5.24) d(x) = a{x)nS . 
Next set 

(5.25) T,(e,)= y f^^iSa), 

dim Sq —j 

for j = 0, 1, . . . , n - A: - 2, where by ^f^l^i and lE23l: 

From the obvious relation 

dim V& = dim Va + dim{V& n V^) 
and the definitions of S and of j (as dim S&), we get dim Va = j + k + 2 , thus 



T,(^,)= U T,,+,.,ASa)ns 



(5.26) =Tj+k+i{ej+k+i)nS . 

Next set 



dimSa=0 
qCck 



and 



— ^^m(eOj • ■ ■ J Em-l) — ( [J '5'a) \ [J Tj(ej 



dimSa=m J<™ 
aCcK 
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for TO = 1, . . . , n — fc — 1. Finally, set 

m<n~k— 1 

From (|5.22|l and H5.2fi|l we obtain the key relation: 

(5.27) n^nns. 

5.4. Proof of Main Lemma. We shall prove this by induction on the dimension 
of S. 

Inductive Hypothesis (n). Let {S,V) be a stratification with dimS* < n. If e — 
(eo, ei, . . . , CdimS-i) decreases sufficiently rapidly and prevents improper tube inter- 
ference, then for x €z S 

clos(x, ft) C a{x) . 

After some preliminary discussion, the proof will consist in proving four claims 
which localize clos(a;, f2) with greater and greater precision. The fourth and final 
claim finishes the induction and hence the proof of the Main Lemma 5.6. 

Repeating here the contents of Definition 15 . 21 for the convenience of the reader, 
the premise of the inductive hypothesis states that 

,/ 2—2 

(1) for m > 0, Cm/em-i < t, with eq < rr, where r = g ^ ; 

(2) for < TO < dimS" - 1 and to G f2,„, if Va C Vi ^ V^, , then \Jse^{^) H 

For a stratification {S, V) induced by the stratification {S, V) we get that tm inherits 
property (1) from em+k+i- Also U8?^(a)) inherits property (2) from Vscrr^+k+ii^) ■ 
i.e., U8?,„(w) n fj^(5a) = , since Vse^+^+^i^) n Te,„_^,^,(55) = whenever C 

Under the assumption that the premise of the inductive hypothesis holds for the 
stratification {S,V) where dim 5 < n, we get from the above inheritance properties 
that the premise holds for all induced stratifications (5, V) with dim S < n — 1. So 
the conclusion of the inductive hypothesis then holds for these lower dimensions. 
This fact will come into play in the final claim. 

To start the induction at n = 1, observe that when dim S ~ 0, S consists of two 
points and the induction hypothesis (1) immediately follows. 

Next assume that the inductive hypothesis (n) holds. To prove the induction 
hypothesis (n+1) let 5 be a sphere dim S = n and e = (eo, ei, . . . , edimS-i) decrease 
sufficiently rapidly and prevent improper tube interference. 

Take x ^ S. If a; e there is nothing to prove. Suppose x ^ ^l. Then 
X S ljj"<n '^ii ^ unique spherical cap Fa{z, ek{x)) which satisfies Lemma 

15.51 Recall Equation (|5.2|l : i.e., 

P«(z,3efe(^)) = (J Fa{y,3ek(x)) . 

(In this expression, as in the convention for spherical caps, the dependence of a and 
.2 on X is indicated by the subscript k{x) on e.) Either z = w G i^k{x), or z ^ ^k(x) 
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and there exists u: € ^k{x) such that diz^io) < Ttk(x) by usmg Lemma l5.8l In either 
case by (|5.6|l we have 

(5.28) P„(z,3efe(,)) cUse^^jM- 

5.4.1. Claim 1: clos(a;,il) C Pa(2, 3efc(a;)). 

Proof: Due to the definition of k{x) and property (2) of Definition 15.21 the sphere 
z,£k{x)) has a non empty intersection with fi: so the distance from x e 
Fa{z, £k{x)) to is less than or equal to the diameter (^^(z, ek{x))) = '^£k{x}- Claim 
1 follows because d{x, S \ Va{z, i^kix))) > '^^k{x)- 

5.4.2. Claim 2: clos(a;, fJ) C 3efc(^)) U (Use,;,, (z) n (7(2)). 

Proof: Let oj S clos(a;,r2). We know by Claim 1 that oj e Pq(z, Se^jj.-))) : so 
we can assume that e Pq(z, 3ej.(2.)) \ \]y,c^^^^{z) n o'(z). We just have to prove 
that oj G Fa{z,'iek(x))- From Claim 1 there exists y S U3ej.(^j(z) fl cr(z) such that 
to G i^a(j/, 3efe(2;)) \ U3e^(^j(z) n ct(z). Project the point u to 3efe(a.)) by using 

the isometry given in Lemma [5. 21 say 

uj' = 0"^(cj) e Fa{z,iek(x)) \ Fa{z,ek(x)) ■ 
Now there are two possibilities. 

Either uj' £ Q, or uj' e Tj,j < k{x) . 
This last possibility does not occur because it would imply 

drelFa{z, <ik(x)) C Tj , 

by using Lemma 15.31 contradicting the definition of k{x) in Lemma 15.51 Hence, 

Lo' G n. 

Assume oj' 7^ uj. Then by the property of geodesic right triangles resulting from 
Definition I5.2r i') we get |w — > jw' — x\ which contradicts the fact that lo G 
clos(a::, Vt). Hence uj = uj' . 

5.4.3. Claim 3: clos(2;,Sl) C dreiFa{z,ek(x))^ (Usefc,,, (z) n (7(2:)). 

Proof: Let uj G clos(a;, fJ). We know by Claim 2 that uj G Fa{z, 'iek(x)^(^3<ik{^) {z)^ 
a{z)) : so we can assume uj G Fa{z, 'itk[x)) \ {z}- Along the great circle on S from 
UJ through X we have a point uj' G dreiFa{z, ek(x)) which is well defined by taking 
the shortest geodesic from uj to uj' . Again their are two possibilities. 

Either uj' £ ot uj' £ Tj, for some j > k{x) . 

The last possibility does not occur for the following reason. Using Lemma 15.51 
we can take j — k{uj') > k{x). This means that uj' G Fa(^^r-^[z{uj'),ek(ui')) for 
some z{uj') G Sa[^,) D Sa with dimS'Q(^/) = k{uj') and dreiF{z{uj'),ek{ui')) n (Tq U 
• • • U Tk{u)')-i) = 0- Also we would have that uj G c/os(w',0), for otherwise uj ^ 
clos{x,fl). Using claim 2, with uj' playing the role of x in that claim, we would 
get UJ G Fa(^^')(z{uj'),3ek(i^'))- Since uj,uj' G S^^^,-^{z{uj')), we would also get x G 
S'^(^,-)(z(w')). Applying Lemma lOl Case HI, we would get x G Fa(i^'){z{uj'), ek{uj')) 
(See Figure IHl) which would imply that k{x) > k{uj'), a contradiction completing 
the proof of claim 3. 
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Figure 8. Schematic Aid to Claim 3. Here the disk corresponds 
to a transversal view of the horizontal strip in Figure [7| 



5.4.4. Claim 4: clos(a:,r2) C a{x). 
Proof: From the previous claim 

dos{x,n) C dreiFa{z,ek{x)) U U3ej^,^)(z) H a{z) . 

li X — z, then, by virtue of (|5.1UI) and the fact that the distance from x to the 
boundary of the efc(a;)-tubular neighborhood of z is greater than ek{x) (see Remark 
15.11) we get, 

clos(a;, n) C Usej^j^j (z) n a{x) . 

So assume x ^ z and 

(5.29) clos(a;, 17) n dreiF^{z, e^x)) \ Use,,., (z) n a{z) ^ . 

Project a: S F^iz, ef^(^^^) along the great circle on 5* from z through x to x G 
dreiFa{z, £k{x))- Since X and x lie on an arc of a great circle which doesn't cross cell 
boundaries, x € a{x). Hence 

(t(x) — <7{x) . 

If a; S ri, then clos(a;, ft) = {x}, and we are finished. If x ^ fJ, we use the inductive 
hypothesis to get 

clos(x, f2) Ca{x). 
We can assert for x e dreiFa{z, e^^)) that 

(5.30) a{x) = a-(x) n dreiFa{z, ek(x)) C aix) 

despite the fact that it does not hold for spherical cap boundaries in general. 
The reason for this is that affine subspaces in V which intersect strata of 5* = 
9reiFa{z, efc(a;)) to form cells in S are precisely those in V which contain V^. From 
(|ST7| and 

clos(a;, rt) = clos(a;, ft) . 

However, we have 

clos(x, ri) = clos(a;, fi) , 
the reason for which is the following. 
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Figure 9. Reduction to a great 2-sphere in three dimensions 
Observe that 

s n s{x, d{x, n))^sn s{x, d{x, h)) 

as illustrated in Figure El where we have drawn the great 2-sphere in three dimen- 
sions determined by the center c and the points z,x,ll! e clos(a;,r2) . S intersects 
this sphere in the meridian containing x and uj. Thus 

cios(x, n)^{sn s{x, d{x, n))) n vl 
= {sc^s{x,d{x,n)))c^^l 

= clos(3;, VL) . 

This finishes the proof claim 4 and also the induction step in the proof of the 
Main Lemma 5.6. □ 

5.5. Conclusion of proof of Proposition l47ll Let S=S(1,0) be the — 1 di- 
mensional unit sphere centered at the origin and V be the set of co-dimensional 1 
subspaces of normal to the edges of the polytope P. 

5.5.1. Proof of Proposition \4-. I}f 3) : the closest points uj^lfltoyQad Ci{vi) are 
in the same cell closure a. Each C(vi) is a polyhedral cone. The half-spaces, of 
which it is the intersection, are determined by co-dimension 1 subspaces orthogonal 
to the edges connecting Vi to neighboring vertices. Let Vp be the collection of co- 
dimension 1 subspaces orthogonal to all the edges emanating from all the vertices 
of P . Let the collection V of co-dimension 1 subspaces of §5 contain Vp. The Main 
Lemma holds for either collection of co-dimension 1 subspaces because the partition 
of S into cells a resulting from V merely refines the one resulting from Vp . That is 
to say that the cell containing clos(a;, fi) given by V is a subset of the cell containing 
clos(a;, fi) given by Vp. This fact is the basis of the proof of Theorem 16. 21 

All that is left of show is that there exists e = (eg, ei, . . . , e„_i) which satisfies 
the premise of the Main Lemma 5.6: namely, e decreases sufficiently rapidly and 
prevents improper tube interference. Upon establishing this, the conclusion of the 
Main Lemma 5.6 will then hold. 
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Let X G S. Because V is a finite collection, the planes in V which do not pass 
through X will have a certain minimal distance to x. Hence, a small enough neigh- 
borhood of X avoids all planes that do not contain x; and from the definition of 
cr{x), we then get that for every x G S there is e{x) > such that the following 
holds: 

(1) for aU y e l],i^^){x) n a{x), Vy=V^; 

(2) if 14 C ^ V„ then U8,(,)(a:) n T,(,)(5„) = . 

By virtue of compactness of the various components of f2fc(eo, . . . , Cfe-i), we can 
then inductively choose a uniform efc for Qk so that: 

,/2— 2 

(1) for fc > 0, ek/ek-i < r, with eo < rr, where r = - — ; 

(2) for every lo £ H Vo^ d V, ^ V^, then Ug., (w) n T,, = , 
as needed. 

5.5.2. Proof of Proposition \^~l][ 2) : for each cell a C C\{vi) there is a > such 
that fl n \Je„{^^) C a. Since d{z, dreiFa{z, ek)) > Cfe, we can choose = Cfe, where 
k is the dimension of a. 

5.5.3. Proof of Proposition \^~l][ 1 ): f2 is not confined to any half-sphere. There are 
several ways of proving this. One way is to chose ek so that the total volume of all 
e-tubular neighborhoods is less than half the measure of S. Another is to use the 
fact that V is the set of co-dimensional 1 subspaces normal to edges of P. Perhaps 
the easiest way is to choose n linearly independent vectors wi , . . . , vn G ^ and 
observe that — Wi, . . . , —v^ are also in Q. □ 

Having completed the proof of Proposition 14.11 we have finished assembling all 
the elements that prove Theorem ll.il 

6. Additional Properties 

Theorem 6.1. The constructed Q is globally absorbing for the family {(/'7I7 € Pq} 
for any Pq a compact subset of the interior of P: i.e., given a compact subset X of 
M.^ there exists an uq (depending on Pq) such that if n > uq then (f)^-^ ■ ■ ■ (f)^^{X) C 
Q /or 71, ... ,7„ e Po- 

We first remark that for a fixed compact subset Pq of the interior of P any 
large enough sphere S centered at the origin will be a globally absorbing region: 
because any larger sphere contracts by a fixed amount under the action of any cj)^ , 
the amount of contraction depending on distance from Pq to dP. However, the size 
of the sphere S varies with Pq. Theorem 16.11 is a stronger result. It states that Q 
can be so chosen that it is universally globally absorbing for all such Po- 

Proof: Let pQ > R where R is given by H3.10|l . We shall show that the set 
Q — PqQoo is globally absorbing. From the strengthening of (|3.5|) to strict in- 
equality which follows from a strengthening of H3.6|) to strict inequality, the other 
participating inequality H3.9I) already being strict, we get that for any 7 G Pq and 
p > Po the set pQoo is mapped by (fi-y into p'Qao where p — p' = 5{p) > is a quan- 
tity which depends on the distance from Pq to dP. The function 5{p) is increasing 
for p > Pq. Thus p' < p — S{po) which leads to 

^71 • • • ^ir^ ipQ) C (p - nS{po))Q C pqQ 
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for some n. □ 

Theorem 6.2. Given a finite set of polytopes and their time dependent dynamical 
systems as defined in Theorem M.lX there exists a Qoo such that the set Q — pQ oo 
is invariant for large enough p for each dynamical system. 

Proof: By virtue of the remark in §5.5.1 we just take V to contain all the co- 
dimension 1 subspaces orthogonal to all the 1 dimensional edges of all the polytopes. 
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